Superfluid response in electron-doped cuprate superconductors 
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We propose a weakly coupled two-band model with d x 2_ y 2 pairing symmetry to account for the 
anomalous temperature dependence of superfluid density p s in electron-doped cuprate supercon- 
ductors. This model gives a unified explanation to the presence of an upward curvature in p s near 
T c and a weak temperature dependence of p s in low temperatures. Our work resolves a discrep- 
ancy in the interpretation of different experimental measurements and suggests that the pairing in 
electron-doped cuprates has predominately d x 2_ y 2 symmetry in the whole doping range. 



Identification of pairing symmetry has been an im- 
portant issue in the investigation of high-Tc supercon- 
ductivity. For hole-doped cuprate superconductors, it is 
commonly accepted that the pairing order parameter has 
d 2 .2_ i) 2-wave symmetry pj. However, for electron-doped 
cuprate superconductors, no consensus has been reached 
on the pairing symmetry. A number of experiments, in- 
cluding the angle resolvedphotoemission (ARPES) 
the Raman spectroscopy and the phase-sensitive mea- 
surements 0, El , suggested that the electron-doped su- 
perconductors also have d x 2_ y 2 -wave symmetry. How- 
ever, the results revealed by other experiments are con- 
troversial 0, H, 0, El El • In particular, the magnetic 
penetration depth data measured by Kokales el al. Q 
and by Prozorov et al. llfj showed that the low 
temperature superfluid density of electron-doped super- 
conductors varies quadratically with temperature in the 
whole range of doping, in agreement with the theoreti- 
cal prediction for a d-wave superconductor with impurity 
scattering. However, the experimental data published by 
Kim et al. [llj suggested that there is a d- to anisotropic 
s-wave transition across the optimal doping. For optimal 
and overdoped samples, they found that the low tem- 
perature superfluid density exhibits an exponential tem- 
perature dependence, in favor of an anisotropic s-wave 
pairing state. 

The above discrepancy indicates that low lying quasi- 
particle excitations in electron-doped cuprates behave 
quite differently than in hole doped ones. To resolve the 
discrepancy, a thorough understanding of the electronic 
structure of electron-doped materials is desired. In this 
regard, the doping evolution of the Fermi surface (FS) 
revealed by the ARPES of m 2 - x Ce x Cu0 4 (NCCO) 
is of great interest. At low doping, a small FS pocket 
first appears around (7r,0), in contrast to the hole dop- 
ing case where the low-lying states are centered around 
(ir/2, 7r/2)|l3j. By further doping, another pocket begins 
to form around (tt/2, 7t/2). The presence of the two sep- 
arate FS pockets may result from the band foldin g eff ect 
induced by the antiferromagnetic correlations |l4l 1 1 5| . It 
may also be a manifestation of the lower and upper Hub- 
bard bands [lf|. At the mean field level, the theoretical 
calculations indicate that these two FS pockets can be 



effectively described as a two-band system This 
two-band scenario is consistent with the conjecture made 
by a number of groups |l7lll8l . ll9| on the existence of two 
kinds of charge carriers in electron-doped materials. 

The interplay between the above mentioned two bands 
can affect significantly the behavior of superconducting 
quasiparticles. A generic feature of a weakly coupled two- 
bands system, as first pointed out by Xiang and Wheat- 
ley |2oj |. is the presence of an upward curvature in the 
temperature dependence of superfluid density p s near T c . 
This intrinsic upward curvature in the superfluid density 
has indeed been observed in electron-doped materials by 
a number of experimental groups 0, 0, E^L Eij ] . Not 
only does it lend further support to the two-band pic- 
ture, but also sheds light on the understanding of various 
controversial experimental observations. 

In this paper, we propose to use a two-band BCS-like 
model with d^ ^2 -wave pairing symmetry to account 
for the low energy electromagnetic response of supercon- 
ducting quasiparticles in electron-doped materials. This 
model, as will be shown later, captures the main fea- 
tures of quasiparticle excitations in the superconducting 
state and gives a unified account for the experimental 
data. Our result suggests that the superconducting pair- 
ing in electron-doped materials is governed by the same 
mechanism as in hole-doped ones, although their phase 
diagrams look asymmetric. 

The two-band model we study is defined by 

H = £ifcct fcg .Cifcg- + Vikk' c\ k ,^c\_ k , yCj^^Cjki 

ika ikk' 

+ (v3fefc'c|fe'T c i~fe'i C2 - fe T C2fc i + h.c)j , (1) 

kk' 

where i = 1,2 represents the band around (ir, 0) and 
that around (tt/2, tt/2), respectively. c\k a and c-ika are 
the corresponding electron operators. V\kk' and V2kk' are 
the reduced pairing potentials for the two bands. V^kk' 
is the interband pair interaction. This model has also 
been used to describe superconducting properties of the 
two-band superconductor MgB2 [2^| . In MgB2 , the inter- 
band coupling is weak since the two relevant bands have 
different parity symmetry |2fij . In the present case, the 
interband coupling is also weak since the strong antifer- 
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romagnetic fluctuations do not couple the first band with 
the second one in electron-doped cuprates. 

In electron-doped materials, the superconductivity oc- 
curs at much higher doping than in hole doped ones. 
However, as shown by the ARPES experiments, the ap- 
pearance of the superconducting phase coincides with the 
appearance of the second band at the Fermi level. This 
reveals a close resemblance between electron- and hole- 
doped materials. It suggests that it is the interaction 
driving the second band to superconduct that leads the 
whole system to superconduct in electron-doped mate- 
rials, and that the pairing potential Vikk' has predomi- 
nantly d x 2_ y 2 symmetry, resembling the hole doped case. 
V\kk' can in principle be different to Vikk' ■ However, 
if pairing in the first band is originated from the same 
mechanism as the second band or induced by the second 
band by the proximity effect, Vikk' should most probably 
have A x i_ y 2 symmetry. 

In the calculations below, we assume that Vikk' (* — 
1,2,3) can all be factorized: Vikk' = 9ilklk> , V 2 kk' = 
92-fk-fk' and V 3 kk' = 9zlklk' , where g%, g 2 , and g 3 are 
the corresponding coupling constants, 7u- = 72/0 = 7fc = 
cos k x — cos k y is the d a; 2_ 1( 2-wave pairing function. Here 
we have implicitly assumed that the first band has the 
same pairing symmetry as the second one. This assump- 
tion can in fact be relaxed. The qualitative conclusion 
draw below does not depend much on the detailed form 
of the pairing function for the first band near (77, 0), pro- 
vided there are no gap nodes on the FS of this band. 

Taking the BCS mean field approximation, the inter- 
action between the two bands is decoupled. It is straight- 
forward to show that the quasiparticle eigenspectrum of 
the i'th band is given by the following expression Eik = 
\J&u + A?7k, where Aj is the gap amplitude of the i'th 
band. They are determined by the following coupled gap 
equations Ai = J2k 7fc(ffi(ci-fexCi feT ) + #3(c2-*4 c 2kt)) 
and A 2 = J2k 7fc(52(c 2 -fcj, c 2fcT) + ff3(ci_jfejcivr))) where 
(• • •} denotes thermal average. 

The above expression of Eik indicates that there are 
gap nodes in the quasiparticle excitations of the second 
band. However, there is a finite excitation gap in the 
first band since the nodal lines of 7^ do not intersect 
with the FS of that band if the system is not heavily 
overdoped. Therefore, as far as thermal excitations are 
concerned, the first band behaves as in a s-wave super- 
conductor, although the pairing is of d x i_ y i symmetry. 
This indicates that the superconducting state of electron- 
doped cuprates is actually a mixture of d-wave and s- 
wave-like pairing states. Apparently, the low tempera- 
ture/energy behavior of quasiparticle excitations is gov- 
erned by the second band since the first band is thermally 
activated. This would naturally explain why the typical 
d-wave behaviors were observed in quite many experi- 
ments 0, IE Q 01 @- However, the presence of the first 
band will change the relative contribution of the second 
band to the superfluid as well as other thermodynamic 



functions. This will suppress, for example, the tempera- 
ture dependence of the normalized superfluid density and 
aggrandize the experimental difficulty in identifying the 
expected power law behavior for a d-wave superconduc- 
tor. 

The superfluid density is inversely proportional to the 
square of the magnetic penetration depth, i.e., p s oc A~ 2 . 
Under the BCS mean-field decomposition, the superfluid 
density of the system is simply a sum of the contribution 
from each band and can be expressed as 

p s (T)=p S;1 (T) + p a , 2 (T), (2) 

where p Si i is the superfluid density of the i'th band. It 
can be evaluated with the formula given in Ref. j2(J. In 
low temperatures, since there is a finite gap in the quasi- 
particle excitations of the first band, p Si ±(T) is expected 
to be given by 

p a , 1 (T)~p i , 1 (0)(l-oe- A i/* Br ), (3) 

where A^ is the minimum value of Ai7^ on the FS of the 
first band and a is a constant. There are gap nodes in 
the second band, therefore p s ,2 should behave similarly 
as in a pure d-wave superconductor and show a linear T 
dependence in low temperatures due to the low energy 
linear density of states: 

P.,2(T)~p. )2 (0)(l-Q. (4) 

Thus, in the limit T <§; T c , the normalized total super- 
fluid density is approximately given by 

Ps(T) _ PsMty T p sA (Q) A'JksT fc -N 

Ps (Q) ~ p s (0) T c p s {Q) K 1 

where p s (Q) = p s ,i(0) + p s ,2(0). 

For a pure d-wave superconductor, as shown by Eq. 
(|4"|) , the slope of the linear T term in the normalized su- 
perfluid density is proportional to 1/T C . However, for 
the coupled two-band system considered here, this lin- 
ear slope is normalized by a factor p s 2(0)/p s (0). The 
zero temperature superfluid density p Si i(0) is a measure 
of the diamagnetic response in the i'th band. It is ap- 
proximately proportional to the ratio between the charge 
carrier concentration and the effective mass in that band, 
i.e., p s ,i{0) oc rii/m*. It is difficult to estimate this ratio 
for each individual band. However, as the FS pocket of 
the first band appears immediately after doping and that 
of the second band appears only after the long range anti- 
ferromagnetic order is completely suppressed, one would 
expect p Si 2(0) to be much smaller than p sA (0). This 
means that p s ,2(0) / Ps(0) <C 1 and the linear T term in 
p s {T) is greatly suppressed. Thus the low temperature 
curve of the normalized superfluid density looks much 
flatter than in a pure d-wave system, although p s {T) is 
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FIG. 1: Fitting curves of Eq. Q to the low temperature 
superfluid density data published in Ref. [Tl| for x = 0.124, 
0.131 and 0.152. The inset shows the doping dependence of 
the fitting parameter A' x . 

still governed by a power law T dependence at sufficiently 
low temperatures. 

In real materials, the low temperature dependence of 
p s (T) I p s {0) will be further suppressed by impurity scat- 
tering and the linear term will be replaced by a T 2 term 
in the limit T < T Q 27] 

P s , 2 (T)^ s , 2 (0)(l-^g-), (6) 

where T is the scattering rate. In this case, p s (T)/p s (0) 
becomes 

Ps(T) ^ Ps , 2 (o) k 2 B T 2 p s ,i(Q) ac - A ' l/fcB T ( ) 
p.(q) ~ Ps (o) 6^r A 2 Ps (o) ■ 1 1 

We believe this formula captures the main feature of low 
temperature superfluid density. Indeed, by fitting the ex- 
perimental data with the above equation, we find that it 
does give a good account for the low temperature super- 
fluid in the whole doping range. This can be seen from 
Fig. where the fitting curves of Eq. J7J) to the mea- 
surement data published in Ref. |ll| are shown for three 
representative doping cases in the under-, optimal and 
over-doping regimes, respectively. 

In electron-doped materials, doping will reduce the dis- 
tance between the FS of the first band and the nodal 
lines of jk- At low doping, the contribution from the 
exponential term is small and the T 2 term is dominant. 
By further doping, A' : begins to drop (the inset of Fig. 

the contribution from the exponential term becomes 
comparable with the T 2 term in certain low tempera- 
ture regime. In this case, the T 2 dependence of p s {T) 
would become difficult to be identified if the exponen- 
tial term is not clearly separated. In heavily overdoped 
regime, the FS of the first band will stride over the nodal 
lines of 7&. In this case, A[ = and p s (T) / p s (0) should 



behave similarly as in a conventional d-wave supercon- 
ductor. This picture for the doping dependence of low 
temperature p s (T) / p s (0) agrees qualitatively with all ex- 
perimental observations. 
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FIG. 2: Illustration of the contributions from the two bands 
to the superfluid density in electron-doped cuprates. 

Close to T c , a positive curvature will appear in p s (T). 
This is a simple but universal property of a weakly cou- 
pled two-band system[2(|. To understand this, let us 
first consider the case g% = 0. In this case, the two bands 
are decoupled and will become superconducting indepen- 
dently. Let us denote their transition temperatures by 
T c °j and T c ° 2 and assume T° t < T c ° 2 . For finite but small 
(73, the superconducting transition will occur at a critical 
temperature close to T c ° 2 , i.e. T c ~ T c ° 2 (Fig. HJ). Just be- 
low T c , p s is mainly contributed from the second band. 
However, when T drops below T^, the intrinsic supercon- 
ducting correlation of the first band will appear in addi- 
tion to the induced one, and the contribution to p s (T) 
from this band will rise rapidly with decreasing tempera- 
ture. Consequently, a clear upturn will show up in p s (T) 
around T^. The appearance of a positive curvature in 
the experimental data of p s (T) , as already mentioned, is 
a strong support to the two-band picture. 

To calculate explicitly the temperature dependence of 
p s in the whole temperature range, one needs to know 
the band dispersion £2^. For this purpose, we adopt 
the expressions first proposed by Kusko et al. |l6( 

&k = ± Ui,k + £i,k+Q ± - £ i,k+Q) 2 + 4 ^ 2 ) /2-^ 

where ± corresponds to the first/second band, Sik — 
—2ti(cosk x + cos k y ) — 4t • cos k x cos k y — At" (cos 2 k x + 
cos 2 k y - 1), t'i = -0.25U and tf( = 0.2U. Q = (tt, tt) 
is the antiferromagnetic wave vector and here 5 is taken 
as a constant, pi is the chemical potential determined 
by the occupation number for each band. It was shown 
that the FS contours determined from this formula agree 
qualitatively with the ARPES data 111! ■ Following the 
suggestion of Ref. 0, 0] , we assume that the second 
band is hole-like. The doping concentration is therefore 
given by the difference x = n e — rih, where n e and 
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FIG. 3: Comparison between theoretical calculations (lines) 
and experimental data (symbols) fl"H for the temperature de- 
pendence of the normalized superfiuid density p 3 (T)/p 3 (0) of 
PCCO at three different doping levels. 

are the carrier concentrations of the first and the second 
bands, respectively. However, it should be emphasized 
that similar results can also be obtained if both bands 
are electron-like. 

Fig. [3] compares the theoretical results of superfiuid 
density for a pure system with the corresponding ex- 
perimental data (symbols) 11] for x = 0.124, 0.131, 
and 0.152. The parameters used are t\ = 5, t% = 
1, (<?i,<72,<?3,n e ,0 = (1.3,1.082,0.005,0.214,0.09) for 
x = 0.124, (1.3,1,0.01,0.231,0.1) for x = 0.131, and 
(1.3,0.984,0.001,0.261,0.11) for x = 0.152. As can be 
seen, the overall agreement between theoretical calcu- 
lations and experimental data is fairly good. It gives 
a strong support to our picture. In low temperatures, 
the theoretical curves exhibit stronger temperature de- 
pendence than the experimental ones, especially for the 
case x = 0.131. This is because the impurity scattering 
was not included in the theoretical calculations. By in- 
cluding the impurity scattering, the linear temperature 
behavior p s will be replaced by a quadratic form. This, 
as demonstrated in Fig. ^ will reduce the difference be- 
tween theoretical calculations and experimental data in 
low temperatures. 

Besides the superfiuid density discussed above, our 
model is also consistent with the phase-sensitive, tunnel- 
ing spectroscopy, and other experiments that support d- 
wave pairing symmetry in electron-doped cuprates. Re- 
cently, the ARPES |l4j as well as the Raman spec- 
troscopy showed that the energy gap is highly anisotropic 
and shows a maximum between the nodal and anti-nodal 
regions. This non-monotonic variation of the energy gap 
from the zone diagonal to the zone axis is not what one 
may expect for a single-band d-wave superconductor, but 
is compatible with our two-band picture. 

More experimental measurements should be done to 



further detect the gap structure in electron-doped ma- 
terials. The scanning tunneling measurement that was 
used for testing the two-band nature of MgB 2 from the 
vortex core state along the c-axis|28|, for example, can 
be used to examine the two-gap picture here. Since the 
interlayer hopping is highly anisotropic [2(| and the c- 
axis tunneling current is contributed mainly from the first 
band, this measurement would allow us to determine the 
coherence length of the first band from the spatial exten- 
sion of the vortex core. Comparing it with the coherence 
length of the second band which can be determined from 
the measurement of H C 2, this will provide a direct test 
for our two-band theory. 

In conclusion, we showed that the temperature de- 
pendence of p s in electron-doped cuprate superconduc- 
tors can be well explained by a weakly coupled two-band 
model. Our work resolves the discrepancy in the inter- 
pretation of different measurement results. It suggests 
that the pairing potential in electron-doped cuprates has 
d x 2_ y 2 symmetry in the whole doping range, same as in 
hole-doped materials. 

We are grateful to T. R. Lemberger, D. H. Lee and 
L. Yu for helpful discussions, and M. S. Kim and T. R. 
Lemberger for kindly providing us the experimental data. 
This work was supported by the National Natural Science 
Foundation of China. 



[1] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000), and references therein. 
[2] T. Sato et al., Science 291, 1517 (2001). 
[3] N. P. Armitage et al, Phys. Rev. Lett. 86, 1126 (2001). 
[4] G. Blumberg et al, Phys. Rev. Lett. 88, 107002 (2002). 
[5] C. C. Tsuei and J. R. Kirtley, Phys. Rev. Lett. 85, 182 

(2000). 

[6] B. Chesca et al, Phys. Rev. Lett. 90, 057004 (2003). 

[7] L. Alff et al, Phys. Rev. Lett. 83, 2644 (1999). 

[8] J. D. Kokales et al, Phys. Rev. Lett. 85, 3696 (2000). 

[9] R. Prozorov et al, Phys. Rev. Lett. 85, 3700 (2000). 
[10] A. Snezhko et al, Phys. Rev. Lett. 92, 157005 (2004). 
[11] M.-S. Kim et al, Phys. Rev. Lett. 91, 087001 (2003). 
[12] N. P. Armitage et al, Phys. Rev. Lett. 88, 257001 (2002). 
[13] T. Yoshida et al., Phys. Rev. Lett. 91, 027001 (2003). 
[14] H. Matsui et al, cond- mat/0410388 (unpublished). 
[15] Q.-S. Yuan et al, Phys. Rev. B 69, 214523 (2004). 
[16] C. Kusko et al, Phys. Rev. B 66, 140513(R) (2002). 
[17] Z. Z. Wang et al., Phys. Rev. B 43, 3020 (1991). 
[18] W. Jiang et al., Phys. Rev. Lett. 73, 1291 (1994). 
[19] P. Fournier et al, Phys. Rev. B 56, 14149 (1997). 
[20] T. Xiang and J. M. Wheatley, Phys. Rev. Lett. 76, 134 
(1996). 

[21] J. A. Skinta et al, Phys. Rev. Lett. 88, 207003 (2002). 
[22] J. A. Skinta et al, Phys. Rev. Lett. 88, 207005 (2002). 
[23] A. V. Pronin et al, Phys. Rev. B 68, 054511 (2003). 
[24] N. Nakai, M. Ichioka, and K. Machida, J. Phys. Soc. Jpn., 
71, 23 (2002). 

[25] I. I. Mazin et al, Phys. Rev. Lett. 89, 107002 (2002). 
[26] T. Xiang, C. Panagopoulos, and J. R. Cooper, Int. J. 



■5 



Mod. Phys. B 12, 1007 (1998). [28] M. R. Eskildsen el al, Phys. Rev. Lett. 89, 187003 

[27] P. J. Hirschfeld and N. Goldenfeld, Phys. Rev. B 48, (2002). 
4219(R) (1993). 



